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Abstract. The Balescu-Lenard equation from plasma physics is widely con- 
sidered to include a highly accurate correction to Landau's fundamental col- 
lision operator. Yet so far it has seen very little mathematical study. We 
perform an extensive linearized analysis of this equation, which includes de- 
termining the asymptotic behavior of the new components of the linearized 
operator and establishing time decay rates for the linearized equation. 
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1. Introduction to the Balescu-Lenard Equation 

The Balescu-Lenard equation is a widely accepted kinetic equation which de- 
scribes the dynamics of a spatially homogeneous plasma. It is 

(1) - V • ( / B(i;, V - v\VF) {F,WF - (VF),F} dv* 

F = F{t,v) is the velocity distribution, F^ = F(t,v*) and (VF)* = V„.F(i, u*). 
(1) is a correction to the spatially homogeneous Landau equation (4). But in stark 
contrast to (4), the kernel B = (Bij) introduces a strong nonlocal nonlinearity: 

V J J\k\<ka \k\^ \e{k,k-v,WF)\^ 
The main advantage of this operator is it's inclusion of the effects of Debye shielding, 
which we will see can lead to very different behavior. Shielding is encoded into the 
collision kernel through the plasma dispersion function: 

(3) e(fc,fc.z;,VF)^l + — -I hm / — f^du. 

In (1), (2), (3), t > and the velocity is w = {vi,V2, fs) G K"^. Also, the wavenumber 
k — {ki,k2,k^) G M^" is cut-off at < fco < co- The parameter, k^, represents 
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the wavenumber beyond which collisions are no longer "grazing collisions" . The 
constant no is the average density, e is the charge, Ae is the Debye length, Ve is the 
thermal speed and rUe is the mass. 

The kinetic equation (1) was derived in a much more complicated form by Bogoli- 
ubov [2,3] in 1946. Then, in 1960, Lenard [10] showed how to write the equation 
explicitly in terms of the distribution, as (1). Lenard argued that this equation 
formally satisfies the expected physical properties including positivity of F, the 
standard conservation laws and the i7-Thcorem. Further, Landau's equation (4) 
was derived as an approximation to (1). Independently and in the same year, 
Balescu [1] also derived (1), which is now commonly known as the Balescu-Lenard 
equation. Due to the nonlocal nonlinearity in the collision kernel, there has always 
been an air of extreme difficulty surrounding this model. As far as we know there 
is only one other mathematically oriented paper on the subject [13] from the early 
1970's. For more information on the physical background and relevance of the 
Balescu-Lenard equation see for instance [1, 5, 9-12, 14, 15, 18, 19]. 

The Landau equation, proposed by Landau in 1936, is one of the most funda- 
mental partial differential equations in plasma physics; see several of the references 
just cited. The spatially homogenous Landau equation takes the same form as (1), 
save that the kernel B(t), v — v*,'VF) is replaced by h{v — v*) = (bij) with 



Above L is a parameter which is logarithmically divergent. L is proportional to 



which is divergent at both zero and infinity and therefore requires a cut-off resulting 
in the well known Coulomb logarithm. The truncation at infinity is needed because 
Landau's equation does not model the efi'ects of wide angle collisions, which is also 
the rational for ko in the Balescu-Lenard case. However the cut-off near zero is 
needed because the Coulomb potential decreases very slowly at large distances. 

The inclusion of the efi'ects of Debye shielding at large distances causes a rapid 
decrease. One of the key advantages of the Balescu-Lenard collision operator over 
the Landau collision operator is that it does not require a cut-off at small wavenum- 
bers and can therefore model very precisely an electrically neutral plasma. On the 
other hand, we will argue that the inclusion of effects at small wave numbers makes 
the difference between these two operators enormous. 

Let's try to make this difference more precise. Roughly speaking, one can see 
the Coulomb logarithm in the Balescu-Lenard kernel as follows 







Switching to spherical coordinates, k = \k\u}, we can write (2) as 




And if L is a constant this is just the Landau kernel (4) because 
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But one of the key results of our analysis shows that the Balescu-Lenard kernel (2) 
can be very far away from the Landau kernel (4) in the following sense. 
Consider the normalized steady state Maxwellian 

^l{v) = (2^)-3/2e-l''lV2_ 

We show that up to some lower order decay, the kernel behaves like 

C e5l''«l'' 

Here vr is v in the direction perpendicular to the relative velocity v — v*: 

(5) K|^=|„p_(^.„)^ 

See Theorem 6 for a precise statement. The main new difficulty in our analysis 

is contained in this observation that the effect of Dcbyc shielding on the Balescu- 
Lenard kernel, when evaluated at maxwellian, is to create an exponentially growing 
velocity factor. Further, since (1) satisfies the H- Theorem, we speculate that this 
exponentially growing factor will be present for solutions to the Balescu-Lenard 
equation (1) at the nonlinear level for large times. We consider Bij{v,v — v* ji) 
above because this turns out to be the kernel of the linearized Balescu-Lenard 
collision operator. 

Before stating our main results, we will linearize (1). For suitable functions F, 
G and H define the Balescu-Lenard collision operator by 

Q[F,G,H]^V ■ \ I B{v,v-v*,\/H){F^VG- {^F)^G}dv* 

We linearize this operator around the normalized steady state Maxwellian. To this 
end, consider the standard perturbation 

F = i^ + ^f. 

Then we can write the Balescu-Lenard equation for the perturbation as 
where the linearized Balescu-Lenard collision operator takes the form 

(6) Lf = {Q[f,^ + Q[^/, ^]} . 

The non-linear part of the Balescu-Lenard collision operator is 

N{f) - {Q[^l + ^if, ^l + ^f, m + v^/] - Vm./, m] - QiVJ^f, m, m]} • 

At first glance, due to the way we have written it, L may seem like a fabricated linear 
operator. But a Taylor expansion of the kernel 'B{v,v — t;*, V{/x -|- s/Jlf}) reveals 
that A^(/) is a nonlinear function of the perturbation /. The terms subtracted off 
on the right cancel with the linear terms in the taylor expansion. 
More precisely, for F satisfying e{k, k ■ v, VF) ^ 0, by (2) we have 

Q[(,,f,,F]=0. 

Therefore, 

Q[f, + ^/Jif,ii + F] = Q[^f, M, F] + F] + Q[^f, ^f, F]. 



4 



R. M. STRAIN 



Let us now look at the collision kernel (2) with (3) for F ^ jj. + ^/Jlf- From (3) 

e{k, k ■ V, V{m + Vm/) = e{k, k-v,\7iJ,) + -^ (^) lim / tlh/E^l^du. 

\k\^ \XiJ vloJm k-[v-u\-ir] 

From here we see the first term in a taylor expansion of B(?;, v — v* , V{/i + ^/Jl,f}) 
in terms of the rightmost term above is B(t;, v — t;*, V/u) and all other terms depend 
on /. Therefore (6) is the linearized Balescu-Lenard collision operator. 
Then the linearized Balescu-Lenard equation is 

(7) i+^/ = 0- 

The presence of physical constants does not create intrinsic mathematical difficul- 
ties. Accordingly, to simplify our presentation, we will normalize all constants to 
one. Because of the null space of L (Lemma 9), a solution to (7) formally satisfies 

/ Mv)[l,V,\vf]^/J4^dv= [ f{t,v)[l,V,\v\^]^/j4^dv. 

We are interested in the asymptotic properties of solutions to (7). 

We will prove time decay to maxwellian in weighted energy spaces. Consider the 
velocity weight 

(8) w{e,d){v) EE (1 + \v\^Y/^ exp + \v\^)i) . 

Above £ e K, (7 > and < < 2. If i? = 2 we further assume < q < 1. 
Depending on our choice of parameters, the velocity in this weight can grow either 
with an arbitrarily low polynomial power or alternatively almost but not quite as 
fast as Define the following weighted norm 

Above di = d,,,. Let (•, •) denote the standard _L^(]R^) inner product. We put a 
zero in the norm to drop the entire weight (8), e.g. ii i = § = then 

\g\l = [ \g\^dv = {g,g). 

In these norms we can show linear decay: 

Theorem 1. Let fo{v) satisfy the conservation of mass, momentum and energy: 

{fo,[^,v,\v\']^ = 0. 
If^ > Q, a solution to (7) with initial data f{0,v) = fo{v) and \ foU < °° satisfies 

(9) \f{t)\o<Ce-^*'\fok, 
where C, A > and p ~ . 

It is our hope that this linear decay will aid in a future in construction of classical 
solutions with small amplitude to the full Balescu-Lenard equation. 

The only previous result for the the Balescu-Lenard equation that we know of 
is the work of Merchant and Liboff [13] from 1973. They show that the spectrum 
of the linearized Balescu-Lenard operator is continuous from zero to minus infinity. 
Additionally they obtain analytic expressions for some spherical harmonic eigen- 
functions. This lack of a spectral gap in the linear operator makes it difficult to 
prove time decay. 
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Indeed, the proof of Theorem 1 requires some development because there are 
many comphcated elements of the Balescu-Lenard operator even at the linear level. 
In Section 2 we will analyse the pointwise behavior of the longitudinal permittivity 
e{k,k ■ v,Vn). We write down another formula for e in (15). Then we use this 
formula to determine the asymptotics of e and thereby show that the Balescu- 
Lenard kernel Bij{v, v — u*, V/x) is well defined (Lemma 3). 

Then, in Section 3, we establish an alternate formula for the Balescu-Lenard 
kernel Bij{v,v — u*, V/x) in Lemma 4 via a series of changes of variables. We use 
this representation to show that the kernel contains an exponentially growing factor 
in Theorem 6. In the proof of Theorem 6 we split the integration region in order 
to squeeze a bit of extra decay out of the kernel, which we later need to use to 
estimate the linear operator. 

In Section 4, we look at the eigenvalues of the so-called "collision frequency", 
which is defined as 

(10) a'^{v) = j Bij{v,v-v*,Wii)ii{v*)dv*. 

Despite the exponential difference in the growth of the kernel as compared to the 
Landau kernel (4), the eigenvalues of the collision frequency decay similarly to 
the Landau case. Using a difi^erent series of changes of variables, we find a useful 
formula for the eigenvalues in and around Lemma 7. Then in Lemma 8 we prove 
that the eigenvalues decay for large |u| as 

\i{v)^ci\og{2+\v\)[l + \v\]-^, \2{v)^C2[l + \v\]-^. 

In the Landau case, the decay is as above if we remove the log(2 + \v\) factor [6]. 
This study motivates the definition of the following weighted Sobolev norm 




If 1? is absent, this means that we drop the entire weight as follows 

\9\l = \9\U - E ^3 ^'^(^) {^^9d,9 + dv. 

These turn out to be anisotropic spaces which are motivated by related spaces used 
by Guo in the Landau case [8]. They help measure the dissipation of (6). 

We study the full linear operator, L, in Section 5. We write down it's null space 
in Lemma 9. Then we split the linearized collision operator as 

{Lf,f) = \f\l-{Kf,f). 
See Lemma 10. The main estimate in our paper is the following 
Theorem 2. For any small r? > 0, there is < C{ri) < oo such that 
\{w'^di{(7'^Vj}gi,g2)\ + \{w'^Kgi,g2)\ < v\9i\tM92\,t,'& + C{r])\gllc(n)\o\g2lc{r|)\o■ 
Here vP' = vP'{(., 1?) and lc7(,,)(u) is the indicator function for the set {\v\ < C{ri)}. 

Due to the exponential growth of the collision kernel (Theorem 6) it was difficult 

to expect that such a result was even true, especially with the possible inclusion of 
these exponentially growing weights. However the exponential growth of the kernel 
is only present in the direction perpendicular to the relative velocity and there is 
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some decay remaing in the other directions. We design a sphtting to show that this 
left over decay is just barely enough to prove Theorem 2. 

In the rest of Section 5, we use Theorem 2 to deduce coercivity of the linear 
operator (6) and exponential decay of solutions to the linearized equation (7). We 
use a standard compactness argument to prove the coercivity of the linear operator. 
But all other arguments in this paper are constructive, including the estimate of 
K. To prove Theorem 1, we use a splitting between velocity and time which was 
previously used to establish decay rates for soft potential kinetic equations in [4, 16]. 

To end this section, we remark that the results in each of the following sections 
build on and crucially make use of the results in previous sections. We begin by 
analysing the longitudinal permittivity. 



The main objective of this section is to determine the pointwise asymptotic 
behavior of the plasma dispersion function (3) evaluated at maxwellian. Since this 
appears to be difficult to accomplish from (3) due to the nonlocal operator our first 
step is to derive an alternate formula. This formula (15) seems to be known but we 
were unable to find a complete derivation in any one reference. Therefore we briefly 
derive (15) for completeness. Then we establish the asymptotics of e in Lemma 3. 

Since we focus on the linear operator, which includes the dispersion function 
only evaluated at Maxwellian e{k, k ■ v, V/u), in the rest of this article we will drop 
the dependence of the maxwellian in our notation, to write only e{k, k-v). Now we 
begin our computation of e at Maxwellian. From (3) we write 



If 7^ 0, let ^^^=111 and {/c^, fc^, fc^} be an orthonormal basis for R^. Define an 
orthogonal matrix Ok such that k ■ OkU = \k\ui as 



2. The Dispersion function, e(|A:|,A; • w, V/x) 




(12) 



Ok = [k^ k^ k^]. 



By this change of coordinates. 




In the last step we integrated out the extra variables U2 and W3. 
To go further, we recall the well known Plemelj formula: 
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which is a distriubtion. Here "P.V." denotes the Cauchy principle value. Thus, 



lim / I , , —dui = lim / — r^du\ 



P.V. / dui+i—e . 



The second term on the r.h.s. is in exactly the form we want. For the first term, 
further decompose ui = — (^-^ — uij + ||| to obtain 

P.V. / du, = ^P.V. / 4 dui - / e-^^/^duu 



We collect the last few computations to obtain e(|fc|,a;) = 1 + |fc|~^\l/(a;), where 
\^ = \bR + i^i and 

= 1 - (27r)-i/2^P.V. / -^dy, 
(13) ^ ^ 



The integral part of '^r is a well studied function in plasma physics [7]. 

We further evaluate via the well known formula [18, p. 184-185, eq. (8.2.13)]: 

dy = ^/2^e-"'/2 / ^y'/^dy. 

-oc x-y ' Jo 

This yields the following simplified version of (13) 

(14) *K(a;) = 1 - a;e-^'/2 P e^'/^dy. 

Jo 

By plugging x=-^ = -^-v = k- v into ^, we can write 

(15) e(|A;|,A;-w) = 1 + |fc|-2*(fc.w). 

Notice that above 5* only depends upon the magnitude of v in the direction of k. 

Now that we have a tractable formula for e evaluated at Maxwellian (15), we 
will study it's pointwise behavior. The key term in the decay of e is ^r: 

Lemma 3. As x ±oo, x^5'j^(x) —1. 

Proof. We first write x^'^r(x) as a fraction: 



x^ 



Next, we apply rHopital's rule to obtain as a; ±oo 



-1. 



^exV2 e-V2 (4-4) 1-4 

That's it. □ 
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Remark 3.1. In particular, Lemma 3 shows that '^r{x), given by (14), is eventu- 
ally negative for large enough x. Prom (15) we have 

Re{e(|fc|,fc-^)} = l+ ^^g,"^ 

Lemma 3 therefore tells us that the real part e will be zero if fc • w is large enough 
and |fc| small enough-|fc| depending on the size of 'i>R{k ■ v). 

There are then lots of places where Re{e} is zero, but if Re{e(|A:|, /c- 1;)} = then 
fortunately 

We conclude that, when evaluated at maxwellian, e 7^ for all finite v. 

By Lemma 3 and Remark 3.1, the kernel of the Balcscii-Lenard collision operator 
(2) evaluated at Maxwellian is well defined. In the next section we will further use 
Lemma 3 to obtain asymptotic estimates of B(t;, v — w*, V/u) in Theorem 6. 

3. The Collision Kernel, B(w, i; - ?;*, V^) 

In this section we consider the Balescu-Lenard collision kernel (2) at Maxwellian. 
Since we will only consider this case, in the rest of the paper we will write 

(16) B{v,v-v*,S/n) = B{v,v-v*). 

We will look progressively more closely at the pointwise behavior of this kernel in a 
series of Lemmas. The asymptotic analysis of the dispersion function from Section 
2 will play an important role. 

We first record some basic properties of the collision kernel which are also shared 
with the Landau kernel (4). A proof of the following can be found in [14]: 

B > and u*B(w, v - v*)u = iS u = c{v - v*) with c e M. 

The main result of this section (in Theorem 6) is our pointwise estimate of the 
asymptotic growth rate the collision kernel, which turns out to be exponential. 
The first step in this direction is to develop a more tractable expression for (16): 

Lemma 4. The Balescu-Lenard collision kernel (16) can be expressed as 

Above, ^jj and are non-negative symmetric matricies that are connected to the 
Landau projection (4) by the formula 

Further, 

(17) E4(-^-<) = E4(-^-<) = o- 

i i 

We define wi{\vr\), w^{\vr\) in (18) and e = i^jj), ^ = i^fj) in (20). 
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The weights are scalar functions which are given by 

wi{M = / smHj{\vR\cose)d0, 

(18) , 

W2{\vr\) = / cos^ 9J{\vR\cos9)d0. 
Jo 

Here J(s) includes the effects of the dispersion function (15) as follows 

do r''° 

(19) Ji^) = 4 / = / Tzm^rrS^TTTi^dp. 



These quantities will help us get all our later estimates. 

The main idea in the proof of Lemma 4 is to use three changes of coor(linatc;s, 
one at a time. The first one is designed to extract the singularity from the delta 
function in (2). The second coordinate change will give us a useful scalar quantity 
in the form of \vr\. And with the final rotation we obtain the integral J which will 
be evaluated precisely in (21) below. 

Proof. Let = v — v* = i^Zy* and {u^, u^,u^} be an orthonormal basis. Define 
R to be the rotation matrix satisfying Rk ■ [v — v*] = k^lv — v*\, e.g. 

R = [ui ^2 ^3]. 



l\k\<ko\k\' \e{\k\,k-v)\^ 



Then we rotate the k variable in (2) with (15), in other words (16), to obtain 
{v,v-v*) = I 

{RkURk)j 5{ks\v-v*\) 
/|fe|<fco 1^1' \e{\k\,{Rk)-v)\^ 
With the formula 5{aki) = ^^{k^) for a > we have 



/ 



\v — w 



|fe|<feo \e{\k\,{Rk)-v)\^ 



We will expand {Rk)i{Rk)j in order to simplify this expression. Note that 



{Rk)i = kiu] + feu^ + feswr- 



3 



Here u = (mj,M2,M3)* with I e {1,2,3}. We therefore have 

3 



{RkURk)j = J2 + E kikmu'yj 

1=1 l^m 

We can therefore write 

Bi,{v,v-vn = ^ '^/g dk. 



1^-^*1 ^|fe|<feo 



\,{Rk)-v)\^ 



Since terms involving k^ vanish because of the delta function, we can define 

^ij = {kiujuj + k2Uiu'j + k\k2 {uluj + . 

This completes the first change of variable. With each new change of variable, 
below, we will elect to redefine as needed instead of repeatedly introducing a 
new temporary notation. 
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Next we evaluate the delta function. Consider k = (fci, /C2)* € and write 

A ki I k2 2 

Evaluating the delta function yields 

i^^-^^i4i<feo \k\'m\,c-v)\ 

where we redefine 



Bijiv,v-v*) = ^ / ^rTTTTl^^dk, 



Before we rotate the k coordinate system again, we look at the following vector 

This is the velocity, v, in the direction perpendicular to the relative velocity, v — v*. 
Since {u^ ,u'^,u^} is an orthonormal basis, we can expand 

\vf = {v ■ u^f + {v ■ + {v ■ u^f. 



Thus, 

This is (5). Also 

C • w = 777 • Wit- 



Now we are set up for another change of variables. 

Let Or be the orthogonal matrix such that {Onk) ■ vr = k2\vR\, e.g. 

Or = [vi vr]. 

Here = {v? ■ v, —u^ ■ w)*. We apply this rotation to obtain 



Bij{v,v- V*) 



_J_ f ^ii 1 

l«-t^l4|<feo \k?\e{\k\,\vR\k2/\k\)\^ -' 



where is redefined as 

^ij = \k\-^{ORk)lulu] + (Oflfc)^«^«| + {0Rk)i{0Rk)2 {u^j + u]u^)}. 
To simplify this expression, we have to write out OrJc. We expand 

{ORk\ = ki {vk)i + k2 {vR)i. 

Hence 

{ORk),{ORk). = kl{v^).{v^). + kUiR)i{VR)j 

+klk2{{VR). {v^). + {VR)j{vji),^}. 

Plugging this into ^ij yields a long expression. But the factor kik2 involves the 
integration of an odd function over an even domain, which is zero. Disregarding 
terms with this factor, we can write 
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where 



(20) 



1 • vf 1 1 {U^ ■ vf 2 2 (m^ • • f) / 1 2 1 2\ 

|„^|2 ^l^J ^ |^^|2 ^1^3^ |^^|2 • 



This completes our reduced expression after a second change of variables. 

For the third and final change of variables, we will split into an angular integral 
and a magnitude integral. The magnitude integral is evaluated precisely in (21). 
For now, we choose polar coordinates as 

fci = psin0, fc2 = pcos9. 

First changing coordinates and second plugging in (15) yields 

_ 1 f^^ {^jjSm^0 + ^fjCOsH}d0dp 

B,,[v,v v)-^^_^^J^ p\e{p,\vn\cose)\^ ' 

1 r2n ,.k„ ^ 1^1^. sjn2 Q + ^2. cos2 9} dOdp 



V-V*\Jo Jo P l + ^+p-4(*| + *f) ' 
1 /"^^ r.l • 2n .2 2/,! Ji\vR\cOSe) 

(4 sm^ e + 4 cos2 0} U0, 



\v — v 



where the last line follows from the definition (19). 

It remains to reduce the integral over [0, 2-jt] to an integral over [0, 7r/2]. Notice 
that '^r{x) is an even function (14) and ^'/(.x) is odd (13), so that ^'j{x) is even. 
In particular, J{x) = J{—x). Since J{x) is an even function (19), the reduction 
follows by first translating ^ w — on the region [n, 27r] and second translating 
6* ^ f - 61 on the region [7r/2, tt]. □ 

Next, it is not hard to evaluate (19) precisely as: 

(21) ^ - fi + + m - f'--'- 



Here we clearly see the Logarithmic divergence of the Balescu-Lenard kernel when 
ko is sent to infinity. With (21), we can determine the asymptotic limits of J: 

Lemma 5. x'^e^-^ ^^J{x) zLVStt as x ±oo. 

Furthermore, if one were to add a cut-off at small wave number in the Balescu- 
Lenard kernel (say fcg > 0) then it would appear in the first arctangent factor in 

(21) as follows. One would replace the factor ^ by ^'^^^ ■ Our proof below 
implies this large exponential growth would then disappear because of cancellation. 

Proof. We will first examine the behavior of the log term. Notice that the only zero 
of ^'/(a;) is at X = and ^'k(O) = 1. From this we see that the argument of the log 
is bomidcd away from zero on any compact set. So there are no finite singularities. 
It is not hard to see that as a; ^ ±oo 

,3,-.V2logA_,^|±^M^Uo. 

^Ux) + ^jix)J 

We therefore only need to look at the term involving the diff'erence of tangents. 
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We now consider the second term of J{x) in (21), which defines the asymptotics. 
By Lemma 3 and (13) 

^i{x) ~ a;2*/(a;) ~ V tt 
Above means the quantities have the same Umit. Similarly 

^i{x) xyn 
Since tan~^(a;) ±577 as a; — > ±00, we conclude 

hm ( tan- (144) - tan- ( '1±14^\ ) = ±.. 



x^±oo \^ \^i{x) J \ */(a;) 

Therefore, 

And the same is true for x^e~^^/'^J{x). □ 

Given u e R^, we define the projection 
(22) PyU = {v ■ u)v 

Then in matrix form we have 

(-^ Pv—v*^ij — ^ij 



V \ V 

\v\ J \v\ 



{v-v*Uv-v*)j 



We now use Lemma 4 and Lemma 5 to get asymptotic bounds for the kernel. 

Theorem 6. Consider the Balescu-Lenard kernel (16) and the relative velocity (5). 
For any < <5 < 1 there exists Cg > such that 

\Bij{v,v- v*)\ < ^^^Cs [1 + \vr\]-'-' . 

Moreover, Vu e and any < g < 1 there exists Cq > such that 

11/ - P pf c 

u''B{v,v-v )u> ^ - l+ t;^ . 

Cg \v-v*\ 

In this sense, the upper bound is exponentially sharp. 

Although the bounds above are sufficient for the rest of our analysis, it seems 
tractable to refine the the bounds established in Theorem 6 and make them sharp 
by evaluating more precisely the integrals in (23) below. 

From Lemma 4, (18) and (20) we see that 

|B,,(t;,t;-OI<4^^^i%^i^^#^. 

\v — v*\ 
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Therefore the upper bound in Theorem 6 requires an upper bound on the weights 
wi, W2- For the lower bound, we use (20) to observe 

tT3f *\ wi{\vr\) (1 W2i\vR\) i 2 

\v — v*\ 

Then using the formula for + in Lemma 4 yields 

U^^U + U'^U = + ^')« = \{I - Pv-v'}u\''. 

We thereby see that for the lower bound in Theorem 6 it is enough to get a lower 
bound on min{wi(|t;fl|), ti;2(|i'ij|)}. This is what we prove. 



Proof. We first establish the upper bound. Prom (18) and Lemma 5, we have 

(23) 



W2{\vr\) < C / C0S2 e -3 d^. 

Jo [1 + I UK I COS 6*] 



61 < f. Then 



Without loss of generality say \vii\ > 1. We split into < < \vii\ ^ and \vn\ ^ < 
3n 

/' 

Jo 



[l + |Ufl|cOS^]3 - \vr\^ [1 + \vr\cOs{\vr\-^)]^ 

We will use 

cosi\vR\-')>l-^\vR\-^'. 
Prom this lower bound and < 5 < 1 we have 

1 + \vr\ cos(|?;fl|-*) > 1 + \vr\ - 1\vr\^~^^ > ^(1 + \vr\). 
Here we have utilized \vr\ > 1. Thus, 

1 e^l^^l' . ^3+, e^'-«l' 



< 2-" 



\vr\' [1 + 1^;^! cos(|t;ij|-^)]3 " [1 + \vr\]^+s ' 

This completes the estimate over < < \vr\~^ . 

Por the second half of the splitting, \vr\~^ < 6 < ^, we have 

/2 .ilDHpcos^e 

d0 < Ce'^\vn\^'=o.^^\^-\-'K 



'\vr\-' [1 + |ufl|cos6']- 
Now we will use the upper bound 

cos(K|-^) < 1 - iK|-2* + ijKi-^^ 

Prom here, we get some weak exponential decay as long as < ^ < 1: 
This is more than enough decay to establish the upper bound. 
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Next we consider a lower bound for min{z/;i(|t)ii|), u'2(|i'ii|)}- By Lemma 5, 

1 r^i"^ gi|t;Rpcos^9 
W\{\vr\) > t; / SVX?Q-r— ■ -3^61 

CJo [1 + |vfl|cos6']'^ 

Jo 



1 



r/2 



W2{\vr\) > jr, / COS^e-— — ■ —d9 

f.7r/2 



Jo 

This time we consider wi and W2 separately. For W2 and < g < 1 we have 

/■7r/2 pcos~^{q) 

/ > / cos2ee5l''«l'^°^'^de> g2cos-i(g)e5l''«l'. 

similarly for wi, 

/■7r/2 /■cos~^(q) / /•cos~^(g) \ 

y >y sin2 6ie5l^«l''^°«'^d6» > ( y sin^ 61^6' J e^l^'^l' > 0. 



These lower bounds for wi and W2 establish the Theorem. □ 

Remark 6.1. It is a basic but important fact that 'B{v,v — v*) = B(u*, w — v*), 
the delta function yields 



Bij{v,v-v*) = I 

= / 

J\k 



kjkj 5{k ■ [v - v*]) ^-j^^ 
\k\<ko \k\^ \e{\klk-v)\^ 



kikj d(k ■ [v - V*]) n / * *N 

^ ^ "-dk = Bij{v*,v- V*). 



l\k\<ko \k\^ \e{\k\,k-v*)\^' 
We will use this later to split up the growth of Bij {v, v — v*) in Theorem 6 between 
\vji\ and Alternatively, the next identity can be used and additionally will be 
useful in other contexts below: 

This is seen by a difference of squares argument: 

p — w*| / \ \u — v*\ J 

In particular, Theorem 6 says that Bij{v,v - v*) e Lf^^{M.l x R^.). 

This completes our estimates for the collision kernel 'B{v,v — v*). In the next 
section we consider the collision frequency. 

4. The Collision Frequency, a{v) 

We recall the Balescu-Lenard collision frequency (10). We will use the reduced 
form of B from Lemma 4 to do an asymptotic analysis of the collision frequency. 
The following quantity will also be used: 

(24) a\v)^{a{v)l}^ = Yl J Bij{v,v - v*f-l^i{v*)dv* . 
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In this section, we compute the eigenvalues Ai(w) and X2{v) of cr(w) in Lemma 7. 
Then in Lemma 8 we transform these eigenvalues into a form which is conducive 
to obtaining precise pointwise information. Somewhat surprisingly, despite the 
exponential growth of the collision kernel (Theorem 6) in comparision to the Landau 
operator (4), the collision frequency still decays at large velocities, and in a way 
which is different but also closely related to the Landau case. 

We also used the collision frequency to define an anisotropic norm (11), | • 
which measures the dissipation of the linear operator L. An equivalent norm will 
be established in Corollary 8.1, using all the asymptotics developed in this section. 
This norm is important for establishing the main results in Section 5. First we look 
at the eigenvalues: 

Lemma 7. The matrix a (v) from (10) has an eigenvalue Xi{v) with eigenvector v 
and a double eigenvalue X2{v) whose eigenspace is perpendicular to v: 

And we can expand a^^v) = ViVjXi{\v\) + {5ij — ViVj) A2(|i'|) with v = j^. 

We remark that the strategy for computing these eigenvalues is to use a series 
of two rotations in a different way from Lemma 4. Here we first evaluate the v* 
integration which is not present in Bij(v,v — v*). 

Proof. First we translate v* ^ v — v* and second we use (12), Ok, to obtain 

^ ' J ' ' JJik\<ko \k\^ \e{\k\,k-v)\^ 

kikj 6{k ■ v*)iJ,{v - V*) 



'|fc|<A 

^-dv*dk, 



\k\<ko \e{\k\,k-v)\^ 

dv*dk. 



kikj 5{k ■ OkV*)^j.{v - OkV*) 



\k\<ko \k\' \e{\k\,k-v)\^ 
kikj d{vl)fj,{v - OkV*) 



dv*dk. 



l\k\<ko \k\' \e{\k\,k-v)\^ 

Next, with k^ = k and {/c^, fc^, fc'^} an orthonormal basis for as in (12), we 
expand the exponent of the Maxwellian as 

It; - OkV*\'^ = \v- v^k^ - vlk^ - vlk^]"^. 

We can further write v = (v ■ k^)k^ + {v ■ k'^)k'^ + {v ■ k^)k^ . Then by orthogonality 

_ Ov*) = (27r)-^/2e"H[('"'=')-''ii'+[(^-'=')-''2l'+[(^-'=')-''3l'}. 

Now we evaluate the delta function and use translation invariance to obtain 

JJ\k\<kAk\' mik-v)\^ 
^ ^ ■ hk, e-i(-'')' 



i|fe|<feo 1^1' \e{\k\,k-v)\^ 
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We will use this formula to compute the eigenvalues. 

Now that we have evaluated the v* integrations via a rotation in the k direction, 
we will simplify the k integrations by rotating in the v direction. Let = v and 
{u^, f ^, u^} be an orthonormal basis for R^. Further define the rotation matrix 



(26) 



= [v^ w^i 



Notice that {Ovk)i = kw] + fevf + k^v^. We rotate the k variable with 0„ in (25) 
to achieve 



{0^k)i{0^k)j 



e 2 



-dk. 



a'\v) = (27r)-i/2 / 

J\k\<ko \ei\k\,0,k-v)\^ 

J\k\ 



'\k\<ko 



m,h\v\/m 



:dk. 



Above, the cross terms in {Ovk)i{Ovk)j disappear because they give you the integral 
of an odd function over an even domain. By symmetry, 



/ 

J\k\<ko 



ko 



zdk 



-j 



:dk. 



f \er-{\k\,kM/\k\)\'- i|fe|<fe„ |e-(|fc|,fci|t;| 
Recall (22). By the spectral theorem I — Pyi = Py2 + P^s, or in another form 



vfv] + vfv^ = Sij-ViVj. 



These last three points yield the result. 

With a sequence of two rotations, we can write these eigenvalues in the form 

2 <5(fc • w*)^(v - V*) 



□ 



JJ\k\<ko ^ ' 

II I 

JJ\k\<ko ^ 



\k\^\e{\k\,k-v)\- 



-dkdv*, 



k V 

W H 



S{k ■ v*)pb{v - V*) 



|fc|2|e(|fc|,fc.^;)|2 

This follows by first rotating the v* variable with (12) to obtain 



dkdv* . 



Xi{v) = (27r) 



X2{v) = (27r) 



-1/2 



-1/2 



(k-vy 



e 2 



\k\<ko |A:|3|e(|fc|,fc-z;)|2 

2 



1 



k\<kQ 



k V 



dk, 



\k\M\klk-v)Y^ 



-dk. 



Next rotate the k variable in the direction of v with (26) to get the eigenvalues 
as written in Lemma 7. These forms of the eigenvalues are similar in form to the 
eigenvalues found for the collision frequency of the Landau equation in [6] . 

We will now transform these eigenvalues into a form which is conducive to as- 
ymptotic analysis in Lemma 8 below. We switch to spherical coordinates via 

A;i=pcos^i, A;2 = psin^i cos^2, fcs = psin^i sin^2, 
< p < fco, < 6ii < TT, < 6^2 < 27r. 
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Then we can write the eigenvalues in Lemma 7 as 



Xi{\v\) = V2^ [ dOx ( ''dp 
Jo Jo 



cos 



2^isin^i e-i^<'="'^^ 



p |e(p, |u|cos^i)|2' 



V 2 7o Jo P |e(p, |t;|cos6»i)|^ 

Plug (19) into the eigenvalues to obtain 



Ai(lt^l) 



dOi cos^ 9i sinfitie-^l^l''^^^' ^'J{\v\ cosfiii), 



A2(|^;|) = £ deism'9ie-i\-\''°'''^J{\v\cose,). 

We change variables as y = cos^i above to obtain 

r+\v\ 



-K 

r+\v\ 
-\v\ 

The fact that the integrand's are even functions yields 



Now wc are ready to look at the decay of these eigenvalues: 
Lemma 8. As \v\ oo 

Moreover, for any multiindex (3 with |/3| < 1 we have 

\DP\i{\v\)\<CHog{2+\v\){l + \v\)-^-\0\, 
\D^\2i\v\)\<C^{2+\v\)-^-\^\, 

where and are positive constants. 

Remeirk 8.1. As \v\ J, these eigenvalues converge to a unique finite limit: 
limAi(a;)=limA2(x) = i,/|j(0). 

xiO xlO 6 \ Z 

J(0) is clearly finite by (19). The eigenvalues thus have no finite singularites. 
Proof. By Lemma 5 and I'Hopital's rule, as \v\ +oo 
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which imphes the decay for Ai in Lemma 8. 
For the decay of A2, we consider 

\v\X2{\v\) = yi^'"' e-'^y'j{y)dy - y^e-^' J{y)dy. 

Since the second term on the r.h.s. converges to zero, as +00 we have 

\v\\2{\v\)^^J\-iy'jiy)dy. 

And the integral on the r.h.s. is finite by Lemma 5. This yields the decay of A2. 

The decay of Ai in (27) follows from taking the derivative of Ai and then applying 
Lemma 5 and I'Hopital's rule, exactly as we have done for the no derivative case. 
The same recipe will establish the decay of the derivative of A2 in (27) . □ 

Then all the results in this section yield the following lower bound: 
Corollary 8.2. 3c > such that 

where = ■w'^{i,'d){v) is defined in (8) and Py is defined in (22). 
Proof. By (22), Lemma 7 and (11), 

\9\l,^ = I w\i,^){M{v){\P,Wg\^ + \v\^\g\^)+X2{v)\[I-P,\Wg\^}dv. 
Plugging in the asymptotics from Lemma 8 yields Corollary 8.2. □ 

An analogous upper bound can also be established in the same way. This com- 
pletes our study of the Balescu-Lenard collision frequency (10). In Section 5 we will 
use these asymptotics to prove bounds for the Linearized Balescu-Lenard Operator. 

5. Compactness of K, Coercivity of L and Exponential Decay 

The main result of this section (Theorem 2) is to show that L can be split 
into as 1/ = —A — K (Lemma 10) where K is "compact" (in the sense of the 
inequality in Theorem 2). This is a standard Theorem for a linearized collision 
operator of a kinetic equation, such as the Boltzmann or the Landau equation. 
However the exponential growth of the kernel (Theorem 6), which is not present 
in the Boltzmann or Landau kernel, creates new difficulties. As a consequence of 
Theorem 2, we deduce coercivity for L in the anisotropic norm (11) in Corollary 2.1. 
Then we finish section 5 by proving exponential decay of solutions to the linearized 
Balescu-Lenard equation. 

First define the projection 

'Pg = {ag + hg-v + cg\v\^]ij}/^{v), 

where ag,Cg € R and hg <E depend on the function g{v). Then we have the 
following standard lemma for the linearized collision operator (6): 

Lemma 9. {Lgi,g2) = (31,^52), L > and {Lgi,gi) = if and only if gi{v) = 
Pgi. In this case LPgi = 0. 
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For the Landau equation, there is a standard argument used to prove this Lemma, 
see for instance [8, Lemma 4]. The proof in the Balescu-Lenard case is exactly the 
same as for the Landau equation because the null space (17) is exactly the same. 

Now we will split the operator L = —A — K using (10) and (24): 

Lemma 10. We split L = —A — K, where 

Ag = di{a^^d,g} + {dia^}g-a^^^^g, 

Kg = -iJ,~'^^'^di^^fi J ^Bij{v,v - v*)iJ,l^'^ i^{djg)* + -^gj^ dv*^ . 

Here and in the proof we use the convention of summing over repeated indicies. 

The proof of Lemma 10 is virtually the same as [8, Lemma 1]. This is expected 
because the Landau kernel and the Balescu-Lenard kernel share the same null space. 

Proof. We define 

Ag = fi-'/^Q[tx,n'/^g,f4, 
Kg = fJ,-'/^Q[^^'/^g,^^,fJ]. 
Then L = -A - K by (6). We will simplify A and K. Notice that 

difj, = -Vijj,, 

dAy^-^l^g\=^-^l^ (d,^'^)g. 

We will use these and the null space (17) of (2) several times below. 
We compute 

Ag = / B,,{v,v-v*){ii,d,{ii^l^g)-{^i^l''g){d,y.)^dv\ 

= ti-'/^di J^^ Bij{v, V - v*)ii.i^'^ [ (dj -^)g + v*g} dv*, 

= pi-^l^dij^^^ B,,{v,v- v*)fi*fi'^^ (^dj + gdv*, 
where we used (17) in the last step. By (17) again, and then by (10) we have 
= /i-'/^di [ Bij{v,v-v*)i^.iJ,'/^(dj + ^)gdv*, 

Next, we take the derivatives on each term to obtain 

= {a'^djg} - a^^ ^d,g + {a^^ ^} d,g + d, {cr^^ | } 5 - | ^5, 

This is the desired expression for A. 



20 



R. M. STRAIN 



Next, for K, we have 



Kg = ^x-^'^^^^^j^^B^j{v,v-v*)[{^l'l^g).^i|l-{^j[|JL^'^g]).|l}dv*y 



By (17) this is 



-M ^^^9i|/Lt j ^Bij{v,v - v*)iiV'^ {{djg)^ + -^g^.) dv* 



And this is the expression we sought for i^. □ 
We are ready to prove Theorem 2. 

Proof of Theorem 2. We first estimate \{w'^{t,'&)Kg\,g2)\. Recall Lemma 10 and 
write 

{w\i,d)Kgug2) 

^{t,'d)Bij{v,v - v*)s/Ti^\{djgi)* + (-ySi) } {dtg2 + ^52} dv*dv, 



w 



+ 



I diw'^{£,^)Bij{v,v-v*),yji^\{djgi)^ + (^gi) ] g2dv*dv, 



where wc recall that (■)* incans that the function in parenthesis is evaluated at 
V*. We will split the integration region several times to obtain the estimate. The 
derivative of the weight function (8) is 

(28) d^{w^i,^)) = w^ii,^)wiVi, 
where 

wiiv) = |2£(1 + Ivfy + q^{l + \vf)i-'^ . 
In particular |wi(w)| < C < 00. Then we can write 

(29) {w^Kgi,g2) = "^2 w^{i,i^)Bij{v,v - v*)y/Ji:ji{d{gi)^diig2dv*dv, 

77 ^K^xIR3 



where 

(30) d{ =d, + ^, dl=d, + ^+w, {v)vi , 

We will estimate {'w'^{£,'d)Kgi, g2) in this condensed form. 

Now we outline the main strategy of the proof. A key point is to get sufficient 
upper bounds for w{£, i?) {v)Bij (v, v — v*)^fi^,fi. We first want to control Bij (v, v — 
v*)y/J^JI by something which approximates the dissipation cr*-' in our norm. But 
we also want velocity decay left over to generate a small constant factor as in 
the statement of Theorem 2. And further we want to show that this bound has 
additional velocity decay to allow us to distribute and control the exponential weight 
(8), which only depends on v. This is done via several splittings. 
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We will now look for a bound for w{£, '&){v) in terms of w{i, ■d){v*). To this end, 
we expand 

1 |2 I |2 f V — V* \^ / V — V* 

= 1^1 " [\v-v*\ ^ [\v-v*\ 

Then, since < i? < 2, we obtain 

el(-M^)* < exp f I (l + |.|^ - .X) 1 eH^- 



If £ > then we similarly have 



2X^/2/ / , ^ 2N 
' — V ' 
r • 



But if £ < then 



(i + N^) ' < i + lt'l^ 



e/2 



Then, by the last few inequalities, under any conditions we have shown 

/ / - * x2\KI/2 . 
w{£, ^){v) < Cw{£, 'd){vR) 1 + ( -v] ) exp ' ^ 



— V \ I I \ 4: \ \V — V 

Since \vii\ = (Remark 6.1) we have 



2\ 1^1/2 

V — V" 

V 



w{1,'&){vr) = w{£,^){v*r) < Cw{£,^){v*) I 1 + . 1^ _ 
The extra factor on the end is needed only if ^ < 0. We have this shown that 

(31) «,(£,^)(.)<^f,^)(.*)(^l+(^^..yj exp(^|(^^..''' 

This estimate allows us to distribute the exponentially growing velocity weight from 

the V variable onto the v* variable. We have to pay with some extra growth of v 
in the direction of the relative velocity, but this can be controlled by terms in the 
upper bound (32) below. 

The next step is to get bounds for \w{£,'d-){v)Bij{v,v — u*)^/i*/i|. Prom Theo- 
rem 6 and Remark 6.1, 

^^exp(-i(^-.)'-i(^-.f) 

(32) \Bij{v,v-v*)./]i:jl\< ^ ^ 



\v-v*\[1 + \vr\]^+^ 

So we have some exponential decay in the direction of the relative velocity, and 
this is how we control the exponentially growing weight (8) and (31). Since either 
< 1? < 2 and < g or I? = 2 and <q <1 there is < < 1 such that 

w{£,'&){v) \Bij{v,v- v*)^iJ,^iJ,\ 

^ Cw{£,&){v*) ( q'( v-v* V l( v-v* 

^-^-^l ^ |^;_„*|[i + |^;^|]3+5«'^Pl 4.\\v-v*\''') a\\v-v*\''' ) ]- 
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We will use this upper bound several times below. 

We first split the integration into a compact region and a large region. For large 
mi > 0, define a smooth function <^i(r) so that (fii{r) = for r < mi and <pi(r) = 1 
for r > 2m\. Then, with (29) and (30), we define Ki so that 

{w'^Kigi,g2) = I yi(|t;| + \v*\)w'^ Bij{v,v - v*)^iJ.<,iJL{d{g{)^(£292dv*dv. 

This is really the hardest term to estimate. We will split the integration region a 
few more times to do it. 

We first split into two regions where (33) yields solid exponential decay in |u| 
and \v*\. Define 



S^ = {\v\>2\v*\}U{\v*\>2\v\}. 



On ^1, if \v\ > 2\v*\ then 



And further 



|v — < + \v*\ < -\v\ 



[v - V*) ■ V ^ \v\^ - \v*\\v\ ^ \v\^ - ^\v\^ 1 



> 

|v — \v — v*\ 

Alternatively, if |?;*| > 2|i;| then 



> 



iv — V 



2 b-t;* - 3' ' - 3' ' 



And similarly 



\v-v*\ < \v\ + \v*\ < -\v*\. 



\v-v*\ ~ \v-v*\ ~ 2\v - v*\ ~ 3' ' ~ 3' ' 



In either case, we plug these estimates into (33) to obtain 



(34) w{i,^){v)Bij{v,v-v*)fx^/^{v)fx^/^{v*) <w{£,'d){v*) 



\v — v*\ 



This is the strongest decay estimate we'll get. 

Here and below we will define terms like {lsiW^Kigi,g2) to be 

{\SiW^Kigi,g2) = / ^ii\v\ + \v*\)w'^ B^j{v,v - v*)y/iJ.^iJ.{d{gi)^dy2dv*dv. 
On ^i we use (34) (and recall (30)) to get 
|(lsiw;'i^i5i,52)| < C 

< Ce / ; r 

J \v — v*\ 



e-lsd^l'+l^*!') 



I'l + jv* \~>m\ 



\v — V*\ 



{w{i,^)digi). \w{i,^)dig2 



{w{e,^)d{gi)^ \w{£,'&)d^g2\ dv*dv. 
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By Cauchy-Schwartz, this is 



e 72 



(kr+k*r) 



1/2 



e 72 



dv*dv 
> 1/2 



\w{£,'d)dig2\ dv*dv 



{w{e,^)d{g^), 



dV 



By (30) and Corollary 8.2, we conclude 

\{ls,w'^Kigi,g2) \ < Ce"™™' |5i|o.,^|52|<t,^. 

Since mi > will be chosen large, this yields Theorem 2 for Ki restricted to Si. 
Next fix < T < 1 and consider the region 



\v — V" 



> \v\ 



U 



\v — V 



> 



Also SI = {\v\ < 2\v*\ < 4|w|}. Then on SI n 52, (33) yields 
(35) w{i,'d){v)Bij{v,v-v*)ij}''^{v)^}l'^{v*) <Cw{e,^){v*) 



\v — v 



where q" = ^-^^ Since r > we can use exactly the same estimates as in the 
previous case to establish 

|(lsfns.w;'(^,i?)i^i5i,52)| < Ce-«"''""|ffiU,^|ff2U,^, 

with q'" = q"'{T) > 0. And this grants Theorem 2 for Ki restricted to Sf n ^2 if 
mi > is large enough. 

It remains to estimate Ki on 5i n 5*2, where 



V — V 



\V — V 



< \vV > n 



V — V 



\V — V 



< \v* 



with < r < 1. Over this region, we do not expect to get anymore general 
exponential decay in and \v*\ out of (32)-(33). However, fortunately. Theorem 
6 allows us the possibility of finding polynomial decay in this region. 

Without loss of generality assume \v\ > 2^/^^^~'^\ which means ^\v\^ — \v\^'^ > 0. 
Then, using (5), (33) and 5*2, with < (5 < 1 we have 



\'w{£,ii}){v)Bij{v,v - v*)y/JIju 



< 



Cw{£,'d){v*) 



1 + \v[ 



{v - V*) -v 



\v — V" 



,-| -(3+<5)/2 



^ Cw{t,d){v*) ^ 2 _ , |2rn -(3+5)/2 



< 



\V — V 

Cw{£,i)){v*) 
\v — v*\ 



-(3+5)/2 
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In particular, since \v\ > 2\v*\ on Sf D S2 we have decay in both variables: 



(36) \w{£, ^){v)B{v, V - v*)^/JI^i\ < 



\V — V 



[l + |^*|2](3+5)/4 • 

And the next step is to use (36) to complete the estimate for Ki on n 5*2. 
By (36) and the definition of t^i we have 

\{lsins^y{t,mm,92)\ 



< c 



f l{|f|>r«l/3} 

isfnsf \v-v*\ 



\w{e,^)4g2\ 



We next use Fubini and Cauchy-Schwartz to obtain 



dv*dv. 



C 



|.|>mi/3 [1 + \v\f+^^''' 
I 



J 

J\v 



If 



< C 



\ 



{w{£,^)d{gi). 



|u|>mi/3 JsfnS" |u - [1 + \v* 



|[l + |^*|](3+5)/2 

2 \ 1/2 

dv* \ dv 



dv* > dv 



l(3+<5)/2^ 



(37) 



We now focus on bounding the term involving {'w{£, '&)d-{gi)t,. 

The main difficulty with this term is controlling the singular factor \v — v*\~^. 
The following is designed to control it. First we remark that on fl S2 

{v - V*) ■ {v - V*) 



\v-v*\ 



\V — V" 



< |t;*r + |u|'" < {l + \/2)\vY. 



Therefore, 5f n 5^ C n {\v - v*\ < 3|t;|^} . Then 

[l + |^*|]-(3+^)/2 



/ 



< 

l\v-v'\<3\v\'^ 

Moreover, 



\v-v*\-^dv* 



\v — v*\ 

V 1/2 / 



{w{e,i})d{gi). 



dv* 



{d{gi)* 



1/2 



Sfnss 



[l+\v*\] 



3+5 



dv* 



l\v-v*\<3\v\ 

Combining these last two calculations yields 



/ \v-v*\-^dv* <C[l + \v\Y. 

two 

i|j;|>mi/3 \jSi 

a 



\v — v 



{w{i,d)d{gi). 



dv* > dv 



ns§ 1 + \v*\ 



1/2 



dv*dv 
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Since v is comparable to v* on S^, this is 

l{|^;*|>mi/6}'U'^(^,'t')('f* 



< c 



[l 



(d{9i)* 



1/2 



nss 



1 + \v* 



3+(5-T 



dv*dv 



And by Fubini's theorem, the integral in parenthesis is 

(digi) " 



I 

J \v* |>mi 



/6 



l + \v 



|3+(5-r 



< C 



\v* |>mi/6 



(rfiffi)* 



l+\v 



*|3+5-r 



|j)-'(;*|<3|?)*|^ 



dt! dv* 



< C 



L 

J \v* I >mi 



/6 



-[l + \v*\f^dv*. 



1 + lw' 



|3+5- 



Since < 5, r < 1 are otherwise arbitrary, we can choose (5 — 4r to be arbitrarily 
close to 1. But this is not needed for the current case, so we merely choose 5 = \ 
and T = I to obtain 8 — At = 4. Then, for mi > large, the integral above is 



< 



C 



my'*log(mi) JvJ^ 



w'^{i,d){v*) 



log(2+|^^1) 



dv*. 



We plug the result of the last few inequalities back into (37) to see that 
\{lsins^y{i,'&)Km,92)\ < 



2 \ 1/2 

dv* 



|f|>mi/3 



w'^{i,'3){v) 



ml^^ ^/\og{ml) 

I 1 2 \ 1/2 

452 



l + \vP+^ 



dv 



< C- 



dv' 



mj^^ log(mi) 



7 «.^^,^)(.)i^^f^l452r<i.)''^ 

By (30) and Corollary 8.2 then 

1/1 2/0 a\T^ \ I , C'l5l|cr,)9|.92|CT,)? 

Klsfnsfw {l,d)K^gi,9i)\ < — nj— -• 

m-{ log(mi) 

This completes the estimate for Ki with mi > large. 

Next define K2 = K — Ki. In this definition of K2 wc have a smooth cutoff 
function < <f2{r) = 1 — <fi{r) < 1 such that (for some mi > large) <f2{r) = 1 
for r < mi and ip2ir) = if r > 2mi. Then from (29) and (30) 

{w'^K2gi,g2) / ip2{\v\ + \v*\)w^B,j{v,v-v*)y/iI^{d{gi)»d2g2dv*dv, 
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where iv^ = {(..'&). Theorem 6 imphes 

= w\£,ir)^2{\v\ + \v*\)Bij{v,V-V*) G L2(r3 X ^3)_ 

Therefore, for any given m2 > 0, we can choose a function ^*) such that 

llCii - <P2'lli=^(R3xR3) < 7 

supp{(/?2''} C {|w*| + |d| < C(m2) < oo}. 
We spht K2 into a "small" part and a "compact" part as 

{K29l,92) = {K2cgi,g2) + {K2s9l,92)- 

Here we have used the splitting 
to define 

(-?^2c6'i,52> = ^ / (P2{v,v*),/Ji:^i{d{gi)»d292dv*dv, 
hi 

{K2s9\-,92) = Y1 j [^'J ~ 'P2]\/J^{di9i)*d292dv*dv. 

We introduce this second smooth cuttofi^ so that we can integrate by parts. We 
estimate each of these terms separately. 

The second term \{K2s9i,g2)\ is bounded by 

llCij - V2 IU^(IR3xM3)||-v/M*M(c^jfl'l)*C^2fi'2||l,2(R3xR3) 



1712 



91 



C 

< — ISlU,^ l52U,i,, 

m2 



where the last line follows from Corollary 8.2. 
After integrations by parts, the first term is 



{K2c9i,92) = / (^^i&{V2{'"^'"*)y/fJ^}) 9iiv*)92iv)dv*dv, 



where 



V- ■ V 

= + y. & = -^vi + + 
Since (w) < C, we therefore have 

I(i^2c5i,52)|<qk2-'llc4/ \9i\'dv\ If \92\'dv\ . 

[J\v\<C{m2) J [J\v\<C{m2) J 

Combining the last two estimates we conclude that for any m2 > 

C 

\{w'^{l,'d)K2gi,g2)\ < — l5'ilo-,,?lfi'2U,^ + C(m2)|5llc(m2)lo|52lc(m2)|o• 
m2 

We thus conclude Theorem 2 for the K2 part by first choosing rni > large and 
then 7712 > is chosen large enough. 

Finally consider 7?)9jC7*gi, 92)\- For m > 0, we split the integration region 

into a compact part and a large part: 

{\v\ < m} U {\v\ > m}. 
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Then by (24) and (27), 

J\v\>m J\v\>m. + PI 



C f Aog{2 + \v\) 



< — w'{i,^) 7 " \g1g2\dv. 

J\v\>m l + 

And then, using Corollary 8.2 and Cauchy-Schwartz, the integral on the r.h.s. is 
< |fl'i|<T,i?|52|<T,tf- Thus, 

f C 

/ w'^{e,i}) \dia^gig2\ dv < -^\gi\aM92\^,i)- 

J\v\>m ™ 

By Remark 8.1, wc see that cr* and dia^ have no finite singularities. Therefore, on 
{\v\ < to}, (24) and (27) imply that there is a constant C{m) such that 

\dia'\ < C{m). 

Therefore, 

/ w\e,^)\dia'gm\dv<C{m)i [ \gif dv] if \gif dv] . 

J\v\<m \J\v\<m J \.^|f|<m / 

We conclude that for any to > 

C 

I (u;2 {£, {>)d,a'gi , 32) | < ^ \9i\„,^ \92\„^^ + C{m.) \gi lc(m)U\g2lc{rn)U ■ 

Since to > can be arbitrarily large, this completes the estimate for the last term 
\{w^{£, i})dia^gi, g2)\ and thereby finishes the proof of Theorem 2. □ 

Next we will use Theorem 2 to deduce coercivity of the linear operator, L, via a 
standard compactness argument. 

Corollary 2.1. 36 > such that {Lg,g) > 6\{I - P}g\l. 

Proof. We use Theorem 2 and the method of contradiction. In this case, we have 
a sequence of functions {gn{v)}n>i satisfying |{I — P}g'„|o. > and 

{Lgr^gn) < ^|{I-P}5„|^ 

Without loss of generality suppose Pgn = and \gn\a- = 1- Consider the inner 
product 

{91,92)0- = y] / cr'^M (digidjg2 + ^^9192) dv. 
Then there exists a go{v) such that 

5n ^ So- 

In other words, {gn,h)cr {go,h)a- as n ^ 00 for all h{v) which are bounded in 
I • lo-. By lower semi continuity, 

I.90U < 1- 

Equivalently \go\l. = lim„^oo(5n) 5o) < lini„^oo \9n\<T\9o\a- = \go\cT- 
Prom Lemma 10 we can write 

(38) {Lgn,gn) = - {dia'gn,gn) - {Kgn,gn)- 
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We claim that 

lim {dia'gn,gn) = {dia'go,go), 

(39) 

lim {Kgn,gn) = {Kgo,go)- 

n— >oo 

These limits will follow directly from Theorem 2. It is then a standard application 

of (39) to prove the cocrcivity. 

We first consider the limit of {dia'gn, g-n)- Splitting 

(dicr'gn, gn) - {dia'go, go) = {dia\gn - go), gn) + {dia'go, {gn - go))- 
From Theorem 2 then 

|(9iO-'ffn,S'n) - (5ia*5o,5o)| < \9n - .9oU (IffnU + l.9oL) 

+C{vi)\{gn - go)lc(r,)|o (|g«lc(r,)|o + |.golc(»?) |o) • 
Then by Corollary 8.2, \go'^c(r,)\o < C \go\^ and \go\^ < = 1. We thus have 

) - {dia''go,go)\ < -^V + C{v)\{gn - 90)^0(71)10- 

Similarly, 

\{Kgn,gn) - {Kgo,go)\ < -^V + C{v)\{9n - 9o)icM\o- 

Furthermore dign are bounded in < C{ri)} because \9n\cr = 1; the Rellich- 

Kondrachov Compactness Theorem thus yields 

lim |(g„ - 5ro)lc(r,)|o = 0. 

njoc 

Then (39) follows by first choosing ij small and second sending n t oo. 
By sending n t oo in (38), using (39), we get 

0=1- {dia''go,go) - {Kgo,go)- 

We therefore conclude 

(Ago,. go) = |.go|^-l<0. 
Since L > we deduce \go\a- = 1- Equivalently, {Lgo,go) = 0. We thus conclude 
that go = Pgo by Lemma 9. On the other hand, since Pgf„ = 0, we also conclude 
that go = 0. This contradicts j^olo- = 1 and thereby establishes the result. □ 

Now that wc have established Cocrcivity of the linearized Balcscu-Lenard oper- 
ator, we are ready to prove Theorem 1. Our arguments in this section are based 
on techniques developed in [16] combined with previous results in this work. 

Proof of Theorem 1. We will prove that a solution to the linearized Balescu-Lenard 
equation satisfying the assumptions of Theorem 1 also satisfies: 

Here we use the notation |/|^ and |/|^ ^ loosely in the sense that we establish the 
above only up to equivalent norms. Then we will show that (40) implies (9). 
We prove (40). Multiply (7) by / and integrate over t; e to obtain 

||/|g + (i/,/) = 0. 
We establish (40) for ^ = ■!? = by Corollary 2.1. 



(40) :^Tl/l^(t) + l/l^,.W<o, 
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Next assume £ & R and < ■& <2. In this case we multiply (7) by w'^{£, 9)f and 
integrate over t; e to obtain 

j^\f\l + {w\l,e)L,f,,f) = Q. 

Our goal is to get a lower bound for d)Lf^ /). By Lemma 10 we can write 

(41) {w\£,'d)Lg,g) = \g\l,, + {d,{w^)a'^ d^g, g) - {w^d,a'g,g} - {w'Kg,g). 
Recall that Ic is the indicator function of the set {\v\ < C}. We claim that 

(42) {w\i,m9],9)>5,\g\l^^-C{ri)\lc^^^g\l, 

where Sq > depends on q from (8) and 77 > chosen small enough. Notice that 
by plugging (42) into the differential equality at the beginning of this paragraph 
and adding the result to (40) for the case 1? = £ = yields (40) in the general case 
(at least up to an equivalent norm). 

To establish (40), it remains to prove (42). We do this first assuming £ gR and 
< 1? < 2. Later we will handle the case 1? = 2 separately. We estimate each of the 
last three terms on the right side of (41). By Lemma 7 and (28) we can write 

di{w^{£,^){v))a'^{v) = w^{v)wi{v)Xi{v)vj. 

Prom the definition of wi{v) in (28), |w^i(i')| < C(l + |i)p)^~^. By Lemma 8, 
\M'")vj\ < Ci2||+pi. Thus for any m' > 

For TO > large, we split the integral into a bounded part and an unbounded part. 

We have 

f 2log(2 + b|) \djg\\g\ ^ ^ ^ ±^ [ ^\og{2 + \v\) ...^ 

J\v\>m ^ + \V? [1 + |^;|]^-^ " J\v\>m ^ + 

<Cm'-^\9\l,,<l\9\%, 

where we have used Cauchy-Schwartz, Corollary 8.2, < 1? < 2 and to > chosen 
large enough. We use Cauchy-Schwartz on the bounded part to obtain 



,log(2+|«|) \d,g\\g\ 



^ TV 2 H , , dv < Cim)\g\^ , 

<l\9\l^ + Civ)\lm9\l 

We estimate the last two terms on the right side of (41) in the same way using 
Theorem 2. This establishes the claim for < < 2. 

For = 2, to prove (41) we will split linear operator in a different way. Specif- 
ically, split L = —A — K (as in Lemma 10) and define M{v) = exp (| (l -|- 
First we can show that there is (5^ > such that 

-{w\£,2)[Ag],g) jj-^ + \v\ya'^v) (di[Mg]dj[Mg] + ^^\Mg\') dv 

-C{6q)\lc{5,)9\l- 
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Second we can establish 

(l + |z;nV^W (d.[Mg]d,[Mg] + |||Mffp) dv > 6,\g\l,^-C{S,) |lc(.,)5|o. 

where 5q = 1 — q"^ — ^ >0 since < <C 1 and > can be chosen arbitrarily 
small. This is enough to establish (42) because the K part is controlled by Theorem 
2. Using slightly different notation, this exact result was shown for the linearized 
Landau collision operator in [16, Lemma 9] in equations (64) and (65) of that paper. 
Since the proof is very much the same in this case, we will not repeat it. 

Now that we have established our claim and thereby the differential inequality 
(40), it remains to show that this implies exponential decay. As in other works on 
time decay problems for soft potentials [4, 16], a key point is to split |/|^(i) into a 
time dependent low velocity part 

E = {l + \v\ <pil+tf'}, 

Then from Corollary 8.2, we have 

where 1b is just the indicator of the set E. Notice that we have ignored the weak 
logarithmic factor. Plugging this into the the differential inequality (40) we obtain 

Thus 



Define A = where for now p = 1 — p' and > is otherwise arbitrary. Then 
(eW|/|2(i)) < Ap(l + t)--ie^(^+*)'|/lBc|^(i). 



pp 

d_ 
di 

The integrated form is 

rt 



|/|g(t)< 6-^(1+*)" |/o|g(i)+Ape-^(^+*)^ / {l + sr-h'^'+^r\fi^^\l^t)ds. 

Jo 

Since |/lBc|g(t) is on £^ = {1 + \v\ > p{l + t)p'} we have 

\flEc\l{s)<Ce-y'^'^^^'''\f\l{s). 
In the last display we have used the region and 

1 < exp (|(1 + l^;^! - |(1 + \v\f^ < exp (|(1 + \v\^)^) e" i''" . 
The integrated form of (40) implies 

(43) \f\lit) < e-^*" (\Mlit) + XplMl j\l + sy-'e'('+^r-iP^ii+sr'ds^ . 

The biggest exponent p that wc can allow with this splitting is p = 'dp'; since also 
p = —p' + 1 we have p' = so that 

1 , ^ 
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Further choose p > large enough so that A = ^ < Therefore the right side 
of (43) is finite. This completes the proof of decay. □ 
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